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Abstract 



In this paper, we consider the Langevin equation from an unusual point of view, that is 
as an archetype for a dissipative system driven out of equilibrium by an external excitation. 
Using path integral method, we compute exactly the probability density function of the power 
(averaged over a time interval of length t) injected (and dissipated) by the random force into 
a Brownian particle driven by a Langevin equation. The resulting distribution, as well as 
the associated large deviation function, display strong asymmetry, whose origin is explained. 
Connections with the so-called "Fluctuation Theorem" are thereafter discussed. Finally, 
considering Langevin equations with a pinning potential, we show that the large deviation 
function associated with the injected power is completely insensitive to the presence of a 
potential. 
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1 Introduction 

Usually, works concerning out of equilibrium stationary systems deal with their local statis- 
tical properties, some assumptions of homogeneity and isotropy being reasonably assumed 
or understood: traditional theory of turbulence, which focuses on local correlations of the 
velocity field furnishes a good illustration of that. A rather new way to study such systems 
was recently proposed by some authors jl], |^, ^| , who preferred concentrate their efforts to 
characterize the process of injection of energy, imperatively required to sustain the stationary 
state. To be more precise, in numbers of such situations, there exists a channel of energy 
injection, usually located at boundaries of the system (rotating blades driving a turbulent 
flow, heated plate in Rayleigh-Bcnard convection, piston shaking granular matter at an edge 
of a vessel, etc . . . ) together with a distinct channel of dissipation, often provided by a bulk 
dissipation mechanism, viscosity or inelastic collisions. This duality is explicitly expressed in 
the dynamical equation for the energy which can always be written in the form E = I — V, 
where T> is proportional to a coefficient of dissipation whereas T is entirely due to the injec- 
tion process: for instance, the evolution of the kinetic energy of an uncomprcssiblc fluid is 
given by 



d t P f v 2 dv 
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where one recognizes easily a surface injection term and a dissipative bulk term (er^ = 
diVj + djVi) . These two distinct "gates" lead to the establishment of a permanent flow of 
energy throughout the system, and is obviously a primordial feature amongst the stationary 
properties of the system. Thus, some experimental measurements jl], ||, || and numerical 
simulations Q were performed to characterize the injection of energy (more easily reachable 
than the dissipation), or more precisely the probability density function (pdf) 7r(e) of e = 
- JI T(t)dt, the averaged injected power during a time interval of length r. In particular, in 
some works |4|, ||, |6| , the quantity 

p(s)= '-log ^ (2) 

was measured, and it was noticed that, according to the conclusions of the so-called "Fluctu- 
ation Theorem" 0, p(e) seemed to tend to a straight line for large r. These observations 
were quite surprising, for the Fluctuation Theorem was established for time-reversal sys- 
tems, a property which is at the heart of the demonstration of the theorem. A convincing 
argument was recently proposed [Q to explain this seemingly universal behaviour: if one 
considers that the signal of I{t) has only finite correlations, large deviation theory predicts 
that log7r(e) ~ t/(e) for large r. Consequently, p(e) = f(e) — /(— e); but for large t, it is 
extremely unprobable to observe large negative occurences of the averaged injected power 
(but not unpossible), since in average, the injected power in a dissipative system is positive. 
As a result, in concrete measurements of /o(e), one can presumably only measure p in a short 
vicinity of zero, for which f(e) — f(—e) ~ 2/'(0)e with an excellent approximation: that is 
why a straight line behaviour is always measured in real or numerical experiments. 

Nevertheless, there was a lack of a real example of a noncquilibrium dissipative model 
in which 7r(e) could be fully computed, in order to show that time-reversal symmetry is 
absolutely required for the Fluctuation Theorem "really" to hold. The initial aim of this 
paper was thus to provide a simple model, where the pdf 7r(e) is exactly computable, and 
which mimics as simply as possible the presence of two distinct channels for energy flow. One 
of the simplest (nontrivial) systems fulfilling these requirements is provided by the Langevin 
equation 

v + jv = ip(t) (3) 
< ip(t)ip(t') > = 2DS(t - t') (4) 

where fluctuations ip(t) and dissipation are considered as two different sources of modifi- 
cation of energy : 



d (I 2 

dt{2 V 



—"/V 2 + ipv (5) 
dissipation injection 



Note that in absence of ip, the system is clearly dissipative and non time-reversal invariant, 
a property also shared by realistic hydrodynamic or granular systems. 

This interpretation of the Langevin equation is clearly uncommon: deriving the Langevin 
equation as an evolution equation for a Brownian particle in a thermalized surroundings, it 
appears clearly that the dissipative term —7?; and the fluctuating term ip are two different 
faces of the action of the reservoir on the particle; the fluctuation-dissipation relation D = 
jksT testifies this profund link. In the present case however, we consider the Langevin 
equation as a given evolution equation, irrespective of its physical origin, and interpret it as 
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a dissipative system it + jv = shaked by a random gaussian force tp; in particular, we will 
make in the following no reference to the fluctuation-dissipation relation just cited. 

Thus, we compute in this paper the pdf of 

e =- f dt'^(t')v(t') (6) 
T Jo 

in the permanent regime using the path integral method (we compute also the pdf of the 
dissipated power). We show as expected that the function p(e) is in this case not a straight 
line (to prevent any confusion, it is worth noting that our result is not contradictory with that 
of Kurchan jjj, who proved the Fluctuation Theorem for the power injected by an external 
operator acting on a Brownian particle: the physical situation considered here is by no means 
the same), though this nonlinear behaviour is extremely difficult to verify numerically. 

But some other interesting and novel features emerge also from our study, which concern 
merely the large deviation function /(e) associated with 7r(e): we show that this function 
is not a regular function but displays a unexpected second order singularity; we discuss its 
physical origin and show that it is intimately associated with the permanent regime which 
allows rare initial fluctuations of the velocity which have deep consequences for the large 
deviation function. 

Finally, we study also the effect of a pinning potential on this large deviation function 
and show that adding a potential does not have rigorously any consequence on the large 
deviation function. This "universality" confirms the relevance of considering the averaged 
injected power as a probe for extracting global features of energy flow into a nonequilibrium 
system. 



2 Free Brownian motion 

2.1 Characteristic functions 

We consider a particle of mass 1, velocity v, whose dynamics is given by (Q), and want to 
compute the pdf of e in the permanent regime. To do that it is convenient to compute first 
its characteristic function 

Sr(A) = <e- A -) (7) 

which is related to the Fourier transform of 7r(e) by FT[ir(ej\(k) = 7?(— ik/r). In some 
works pj, one computes already at this stage the asymptotic exponential dependence of 
the characteristic function: (e~ AT£ ) ~ e T9 ' A ^, and retrieves /(e) as the inverse Legcndre 
transform of g(X): /(e) = g(X) + Ae, </(A) = — e. We will see that this procedure is not 
appropriate here, for reasons which will be made clear later. One prefers thus compute first 
exactly 7?(A), what is here fortunately feasible. 

Let us consider the stochastic equation (ft). The propagator P(vi, t\vq, 0) can be easily 
expressed in terms of a path integral jll], [T^rp^] : 

P( Vl ,r\v ,0) = d T x [ { \w}exp(-^-[ dt(v + jv) 2 ) (8) 

Jv(0)=v V ^ U JO / 

From this formula, one can deduce the probability density associated with a given path: 



V([v(u), < u < r]) = exp - -L J dt (v 



iv? (9) 
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As a result, one has 

f v(t)=v 



dv x \ [Vv]V{[v]) exp -A / v(v + jv) (10) 

-oo Jv(G)=v \ JO / 

/oo pv(t)—vi / i /*r > 

d Wie -(l+A)("?-^) / [x>«] exp -— / [v 2 + ( 7 2 + 4DX-f)v 2 } 
-oo Jv(0)=v o \ ^JO 



(11) 

where (. . . }„ designates an average over the realizations of v such that v(0) = vq. The path 
integral in (|ll|) is well-known and its value can be exactly computed pi[: 



2 i 2 2l . 

orv \ = 



/ [^w]exp - — / [* 



^ sinhar) exp f-EL ^ + <*) cosh ^) ~ ^ ) (12) 
/ P V 4£> sinhar / v ; 



Thus, let us define 



7 = 7 r (13) 
A = 2L»A/ 7 (14) 



7 ? = V / 1 + 2A (15) 



one has 



1/2 

7r W0 (A) = e^ 2 (cosh 777+ — sinh777 J x exp f ^ ^ J] (16) 

\ V J \2D r] coth 777 + 1 + A / 

To get 7?(A), one has to average over vo with a probability oc exp(— Vqj/2D), since it is the 
distribution of v(t) in the permanent regime: 



(~ ~ \ -1/2 
1 + A — A 2 /2 \ 
cosh 777 H sinh777| (18) 
V J 

The calculation leading to (|l6|) and (|l^) assumes a priori A > 0, but it is useful to study 
the properties of the analytical continuations of these formulae. First, one remarks that the 
square root in the definition of 77 does not induce any breaking of analyticity, for it appears 
always in quantities such cosIi(t77) or sinh(?77)/?7 which are entire functions of A. Thus, one 
will assume in the following that 77 = x + iy has positive real part. 

Let us first look at 7r„ (A) for vq = (the exponential term will not modify the analytical 
properties of tt Vo (A)). It is astute to write it as[] 

(e- XT % = (i±^) x (l + \ (r> + i) tanhr^) (19) 



this manipulation makes the leading term in the limit of large r explicit and allows for a simple localization 
of the cuts. 
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and assume again a cut in the negative real semi-axis for the definition of the square root. 
Breaks in analyticity can arise if only one of the two terms of (|l9|) (separated by "x") is not 
defined (a superposition of two cuts restores the analyticity, since they are associated with 
a square root). The last term is not defined for 77 = x + iy such that 

x = and 1 < ^ ( y - - ) tan jy (20) 



y 



This induces in the A-space a dashed half cut localized in the negative real axis, beginning 
at the value A_ of A less but closest to —5 such that 




2A_| ; = ) tan7V|l + 2A-| (21) 

1 + 2AJ 



As 7 oc t, it is clear that A_ — > — \ when r — > 00. 

The first term e^ 1_?? ^/ 2 \j2j (1 + e~ 2r w) , as a function of 77, is analytical in the region 
x (except at points x = 0, y = f-[7r]). But if one considers it now as a function of A, it 
appears cuts in the A-space due to the prescription x = Re(7y) > 0. It is easy to verify that 
these cuts arc defined by 

(1 + 2A) G i[-(37r/2 + 2fc7r) 2 ,-(7r/2 + 2fc7r) 2 ], for k G N (22) 

7 

and again are located in the A < — ^ half axis. 

To summarize, n Vo W continued on the whole complex plane (as a function of A) has a 
cut, dashed line shaped^ and localized on the negative real axis. It begins at a value A_ less 
than —5, but tends to — ^ for large r. 

For 7?(A), the situation is similar, and gives also this dashed negative cut. But there is a 
fundamental discrepancy, for the "second term" [1 + (1 + A — A 2 /2) tanh(777)/27y] -1 ' 2 gives 
here a novel cut, plain and localized on the positive real axis, and beginning at a value A_|_ 
solution of 



1 + 2A-, 



-6\/l + 2A + I x tanhTV 1 + 2A + = 8 (23) 

; 1 + 2A +/ 



When t — > 00, one has simply A+ — > 4. One has summarized these analytical properties on 
figure |l|; it is worth noticing that the extra cut of 5r (A) has deep consequences on the shape 
of the large deviation function, as we show in the following. 



2.2 Pdf of injected power and large deviation functions 

Why did we study analytical properties of tt Vo (A) and 5r (A) ? A priori, we could have remarked 
that 7r„ (A), as well as tt(A) (we will use henceforth the notation 7r(„ )(A) to designate both 
7f„ (A) and 7r(A)) are such that 

7? M (A) ~ exp[r 5 (A)] (24) 

T — ► OO 

with g(X) = l(l- \Jl + 2A j (25) 

2 We do not give further details on the precise structure of this "hacked" cut, for they are of no importance in 
the following. 
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a) 




i Im(A) 






Re(A) 



b) 



A Im(A) 



Re (A) 



Figure 1: Cuts of functions (a): tt V0 (A) and (b): vf(A) in the complex A plane (see text for 
details). 

(where "~" means an equivalence between logarithms |To| ). and using a traditional recipe, 
obtain the large deviation function of 7r(„ ) (e) as the inverse Legendre transform of <?(A). This 
is actually correct for tv Vo (e), but gives wrong results for 7r(e), because the above mentioned 
procedure neglects completely the presence of the extra cut in n(e), whose origin is the 
prefactor of the exponential leading term. 

Thus, it is more suited to first express the Fourier inversion of 7f(„ )(A) properly, and only 
thereafter extract the associated large deviation function from a saddle point expansion JTij . 

Let us define e — e/D the dimcnsionlcss injected power. One has 

7r (oo )(e) =D xtt( Vo) (b) (26) 
= H dX %o)(7A/2 J D)exp^eAj (27) 

The leading exponential term in this integral is exp[rh(X)] with h(X) = ^(eA + 1 — V 1 + 2A). 
The saddle point expansion method requires to distort the usual path of integration A G iM. 
in such a way that h(X) be always real. Let us parametrize yl + 2A = x + iy (with x ^ 0). 
Two paths ensure Imh(X) = 0; they are 

M x ) = x 7 1 ' ( x G I - 00 !) and ( 28 ) 



2 

A 2 (y) = ^f^-y 2 -l)+^ , (ye]-oo,+oo[) (29) 



1/1 „2 ,\ , :V 



E 



(see figure H) and give /i(Ai(x)) = ^[|(x 2 - 1) + 1 - x] and h{X 2 (y)) = ~%[§y 2 + ^(e - l) 21 
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Figure 2: Paths Ai and A2 where equation Im[<?(A)] = is verified for s > 0. One indicates also 
the schematic behaviour of Ke[g] along these paths : if it grows to 00 ("up") or decreases to —00 
("down"). The saddle is located at the crossing of the two paths. Note that A2 does no longer 
exist if e ^ 0, but in that case Ai vanishes at infinity. 



respectively. 

Owing to the prescription ReVl + 2A 0, the path A2 does exist only if e > 0. The 
next step consists in choosing the right path of integration. From now, the computations for 
■K Va (e) and ir(e) differ. 

For n Vo (i), the situation is relatively simple, since neither Ai nor A2 crosses the cut-off 
of 7? „ . In the following, we give details only for Vq = because it simplifies a bit the 
computation; we postpone remarks concerning the incidence of a non-zero initial velocity. 

It is easy to see that A2 is a valid integration path if e > (in particular, the prefactors of 
the exponential do not cause problems of convergence). If e < 0, the path A2 is not defined, 
but using a semicircular contour, one shows easily that (evidently) 7To(e) = strictly in this 
case. To summarize, and after some calculations, the pdf 7To(e) of the dimensionless injected 
power e = dtip(t)v(t), knowing that the initial velocity vo is zero, is given by 

w [ if e< ( ' 

with 



1(e) = x / dri 77 e 6 ^" ( - I X I 1 H — I n-\ — ) tanhr?7 

4i7r J a _ ioo \ 2 ) \ 2 V V 1 



1 / 2 / 1 / 1 \ \ -V2 



l + e- 2 "^\ 1 1 / 1 



(31) 



(cr is any positive number). It is quite difficult to simplify this prefactor. In the large r limit, 
it is equivalent to 

1(e) ~ i— , (32) 

v V^ooV ire{e + 1) v ; 

but this equivalence is not uniformly valid for all e. In particular, for fixed t, the large i 
values give a e -3 / 2 regression instead (a regime which arises for e > 7). 
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From (pfl) one gets the following large deviation function 



/„( £ >-i( £ --l) 2 xfi( £ -) (33) 

(9 is the Heaviside function), what matches the rapid evaluation above mentioned. 

How are modified these results if vq ^ ? Essentially, an "energetic" initial condition 
gives rise to a small interval of possible negative power injection. To be precise, if Vq =/= 0, 
the probability 7r„ (e) is no longer zero if e s] — |uo|/2r, 0[ (of course the positive part of the 
pdf is also slightly modified). But these modifications are minor, and in particular, they are 
unable to affect the associated large deviation function (the negative window vanishes when 

T — > OO). 

Let us now look at 7r(e). As the function h is the same for both tt Vo (s) and 7r(e), the 
paths Ai and A2 are similarly defined. The essential difference comes from the positive cut in 
the analyticity of 7?(A), which can cross the steepest descent paths Ai and A2 (see figure 
If e > 1/3, the path A2 avoids the cuts and can be chosen as a valid integration path (fig 
(a)). On the contrary, if < e < 1/3, it crosses the positive real cut and a portion of the cut 
must be crawled along to close the path (see fig. || (b)). If e < 0, the parabola A2 does no 
longer exist, but Ai is valid (more precisely a U-shaped path sticked on each side of the cut 
- see fig. H (c)) and leads to a non zero result for the probability. After some computations, 
one can deduce the following result: 

J(e) x e-i^" 1 ) 2 if e > 1/3 

*r(e) = { J(e) x e"*^" 1 ) 2 + K(e) x e^-i) if < e < 1/3 (34) 
K(e) x e^ 2i ~V if e < 



i rtr_ j-Kn-^r ~ifi_±_ ,„< _ ^ 



with 



- /-1/Max(0,e) _ / 1 , -2m\ 3 / i 



(35) 



+ 



K(i) = e—r^ £ - 2 >^- dx x e^ x ~ x i I — ^ J x I — (x 4 - 6^ - 3) tanhx7 - 1 



(36) 



where x+ is positive, defined by A+ = (x\ — l)/2 (note x + — > 3 for large r). For large values 
of t, one can simplify a bit these formulae: 



T — too ITT 



er-Hoo 



a— too 



m ~ - dnne^-^J - J ■ (37) 



(3-^(77 + 1)3 



Ri£) r~°oVj Q ^ U + lJ ^ (38) 

but, as for /, the latter simplification for J is not uniformly valid for all concerned values of 
e. 

It is easy to verify that neither J nor K have an exponential leading term; thus, the large 
deviation function is easily computed as 

m = l -S (f - 1)2 if ^ 1/3 (39) 

1 7(2e-l) if e < 1/3 



8 



(a) e>\ 





N 'ir 








Saddle 









(b) 0<£<± 





















-i 

i 


Saddle 



(c) e€0 



Saddle 



Figure 3: Sketch of different real paths, according to the values of e. The positive real cut fixes 
the location of the saddle at his extremity A + as soon as e < 1/3. 
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Figure 4: Semilog plots of the pdf of injected power ir(i) for several values of 7 = 7T. One plots 
together formula ( |34| ) (line) and results of numerical simulation (dots). 

This bipartite shape of the large deviation function is closely related to the location of the 
saddle of the integration path, which remains fixed in the complex plane as soon as e < 1/3 
(see figure |J). These results are successfully compared with numerical simulations, see fig. 

I 

2.3 Discussion 

These results address naturally some questions, in particular concerning the surprising struc- 
ture of /(e): this function displays a second order singularity located at a odd s = 1/3 value. 

Let us first give some remarks on the shape of 7r(e) (we restrict the discussion henceforth 
to situations where r 3> 7" 1 )- This pdf is a rather asymmetric curve, what can appear at 
first sight surprising : as the renewal of the noise ip is independent of the particle veloc- 
ity v, occurrences of positive or negative instantaneous power injection ipv arc completely 
cquiprobablc. Actually, the long time interval during which the mean is performed is of 
crucial importance and is responsible for the peculiar shape of the pdf; to understand this, 
let us consider an occurrence of a (rare) large positive fluctuation of injected power: this 
occurrence understands that a favourable sampling of the noise is realized, so that very often 
the noise gives energy to the particle. Consequently, during the process, the energy has a 
global tendency to increase, as well as typical values of the velocity (despite the always acting 
dissipation); as the injected power is directly proportional to the velocity, one sees that the 
direct effect of the positive injection of energy is to enhance typical values of v implied in 
the evaluation of e: this favourable feedback makes finally the occurrence of the considered 
fluctuation more likely, since less efficiency of the noise is globally required to generate the 
fluctuation. 
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On the contrary, let us consider a (large) negative fluctuation: the scenario is here in- 
verted, since the typical velocity will certainly decrease during the mean process : to ensure 
the large expected value of power ceded back to the bath, one is then compelled to begin the 
motion with a large kinetic energy, what is exponentially improbable : the feedback is here 
clearly unfavourable, and diminishes comparatively the frequency of such occurrences. 

This primordial role of the initial energy on the negative power injection mechanism gives 
actually also the explanation for the presence of the negative tail in the large deviation func- 
tion of 7r(e): each fixed initial velocity is surely relaxed within a characteristic time oc 7" 1 , 
but this time is longer the higher initial velocity Vo; thus, when this velocity is statistically 
distributed according a Gaussian, rare large initial velocities construct the negative tail of the 
distribution as well as the associated large deviation function, which precisely characterizes 
rare events. It is worth to note that the specific negative tail due to the thermalization of vq 
begins actually at the value 1/3, and thus affects substantially also the positive part of the 
distribution. 

This singularity of the large deviation function can be interpreted as a phase transition, 
if one writes tt (A) as a functional integral over the realizations of the noise tJj. Integration on 
v gives 

7?(A)oc /o^expf-i / ip 2 -^- f dt dt' ^MOe" 7 ''"*' 1 + < ^^~ \ / dt ^(t)e" 7t 

(40) 



4 Jo 2 J 2 7 



and this expression can be viewed as a configurational partition function of an unidimen- 
sional line ip(t) of length r confined in a quadratic potential V[ip] (x f ip 2 , with short range 
homogenous interactions (second term) and an additional local destabilization term (third 
term), which is a direct consequence of the thermalization of vq. If A is positive and too 
large, the combined effects of the short range interaction term and the third term -which 
can be viewed more or less as an inverted parabolic potential acting in the vicinity of the 
t — end of the chain only- destabilizes the chain which is no longer confined by V. 

We would like to make here a little mathematical digression. One can compute the char- 
acteristic function from the preceding formula. This leads to the following formal expression 




1 ^ / 2A 
- > log 1 + 

2 ^ & \ 1 + : 



2A 



,•2 



--logf 1-4A 2 V — — : ^ = I 

2 *\ ^ 1+^2 + 7(1 + ^) 1 + 2A + 4/V 

where x n is solution of the implicit equation jx n + 2Atanx n = mr. That this involved 
expression is exactly the simple formula ( |l§| ) is quite amazing, and one can derive some 
nontrivial mathematical relations from this connection; to mention but a few, one has for 



instance 



a relation which is easily checked for 7 — > {x\ ~ \[2pf) and 7 —> 00 (the sum becomes a 
Riemanian sum). 
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Figure 5: Injected and dissipated power pdf for the free Brownian motion (p/r = 10, D = 1). 
The two distributions are very similar but in the vicinity of zero, partly due to the constraint 
"yv 2 0, whereas ipv ^ 0. 



To conclude this paragraph, it is important to note that, in contradiction with intuition, 
initial conditions can have deep consequences on the shape of large deviations functions, 
even if the process displays finite time correlations, provided that these initial conditions are 
statistically distributed over an unbounded interval, what is almost always the case if for 
instance stationary processes are considered. 



2.4 Dissipated Power 

The pdf of dissipated power can also be computed along the same line of reasoning. One 
gets for the characteristic function of the stationary process: 

^exp— A7 J dt v 2 (t)^ = e?l 2 (cosh rfy + T/sinhr/7) 5 (43) 

and one sees that there is not any positive cut in that case. As a result, the large deviation 
function of the stationary dissipated power 5 = ^ fj" jv' 2 is easily derived as fdissip{5) = 
fo(S). Consequently, one expects the pdf of dimensionless injected and dissipated power to 
be (in the large r limit) similar but in the zero injection (or dissipation) region. This is 
effectively the shown on figure |l 



2.5 Note on the Fluctuation Theorem 

As explained in the introduction, our model is a good system to test the possible universality 
of the conclusions of the (Evans-Cohen-Morris) fluctuation theorem, since the exact result 
is at hand. From (|3^), one has 

1 n(e) { 4 7£ if e < 1/3 

,® = -]og-2E^ ~ 7 3 7 (44) 

T 7T(-£)t^oo --ye + -<y — if £ > 1/3 
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This function is clearly not a straight line, as it would be if the FT held here (not even the 
slope at e — is in accordance with the FT. which would predict p{e) = je). Thus wc 
exhibit here an example where the conclusions of the theorem are not verified (due to the 
simple fact that the situation considered here does not fulfill the hypotheses required for 
the Fluctuation Theorem to hold). On this "negative" result can we make two comments: 
first, it is not contradictory with those of Kurchan One looks for the power injected by 
a random fluctuation in a dissipative system, whereas the Fluctuation Theorem established 
by Kurchan considers the power injected by an external operator in a system in equilibrium 
with a thermostat. We already noticed that the first situation is much more appropriate 
to describe realistic systems driven far from equilibrium. Second, formula (^4|) illustrates 
well the fact that Fluctuation Theorem seems to hold in so large a number of experimental 
situations, as explained in Q|: in the vicinity of e = 0, p(e) must always have a straight line 
behaviour, as a consequence of the large deviation law; on the other hand, as large negative 
values of e are extremely unprobable when t is large, it becomes practically unpossible even 
to only measure p{e) for large e and large r with enough statistical resolution: possible 
deviations from the straight line are just even not measurable. In our case, crossover occurs 
for e = 1/3 and for this value, 7r cx exp(— 57r/3) which is of order 10 -8 only if jt = 10 . . . 
Our model is thus a good illustration in favour of arguments given in [|] against an universal 
applicability of conclusions of the Fluctuation Theorem. 




0.1 0.2 0.3 0.4 0.5 0.0 0.7 0.8 0.9 1 

e 



Figure 6: Function r _1 \og[Tr(e) /ir{—e)] for jt = 8. Circles come from numerics with two millions 
points of statistics, and stop before the crossover of the theoretical curve (solid) due to the lack 
of negative points with large absolute value. 



3 Confined Brownian Motion 

The reasoning concerning the asymmetry of ir(e) suggests that certain characteristics of this 
pdf seem — in the limit of large r only — to be independent of the details of the particle 
dynamics, since it is based only on considerations on energy and its conservative character. 
These observations led us to infer a possible insensivity of tt with respect to other micro- 
scopic times than 7 _1 (7 _1 itself cannot be neglected, since this time plays a role in the 
process of dissipation of energy; and indeed, the curves for different 7 have different and non 
superposable shapes), in cases where the initial system would have been complexified. Thus, 
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we considered several confined Langevin systems: 



x + ~fx + V'{x) = V>(*) (45) 
< ip(t)il>(t') > = 2D5(t - t') (46) 

where we used for V{x) an harmonic potential V{x) = ^uj 2 x 2 , a non linear "hard" potential 
V(x) — ^uj 2 x 2 + jx 4 (a > 0), a non linear asymmetric "soft" potential V(x) = uj 2 (a 2 e~ x / a + 
ax), and also a bistable (p 4 potential V(x) — Vb(x 2 — l) 2 . We numerically calculated for all 
these potentials (as well as again the free case for comparison) the pdf 7r(e), for large values 
of r (7 = 1 for convenience). 




e 

Figure 7: Injected power pdf for several trapped Langevin particles with D = 1, r = 20 and 
7 = 1: three harmonic potentials (free motion u = 0, u = 1 and u = 10), hard and soft 
nonlinear potentials (u> = l,a = 3 in both cases), bistable (^-potential with Vf, = 1 (strong 
anharmonicity is explored in this case - see details in text). 

The results (see fig. ^ are extraordinarily surprising, since one cannot distinguish the dif- 
ferent curves from each other! One must keep in mind the fact that the dynamical behaviour 
of x(t) is completely different in all these cases: fully isochronic or anisochronic oscillatory, 
overdamped, bistable, diffusive, these different Brownian motions lead all to apparently the 
same curve, a coincidence which goes beyond all expectations. 

To understand this phenomenon, let us look at the characteristic function 7r ro (A), where 
i"o = (xo,vq) designates the initial conditions. Of course, one cannot compute it exactly, 
since in general the dynamics is nonlinear. Nevertheless, it is possible to express it in a 
fruitful form. 

From p3| , one can derive the path integral representation of the solution of the Kramers 
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equation: 

r ( r )= r i 



P 7 (r 1 ,r|r o ,0) = e^ 2 x f [Vx] exp (--L [ dt [x + 7 i + V'(x)]A (47) 

Jr(0)=r \ 4L» J / 



'r(0)=r o 

(the index 7 recalls the value of the damping). As before, one gets from this 

r( T )=n 

[Z 

'r(0)=r o 



7r ro (A) = / dn / [2?i] exp (l--L [ dt [x + jx + V'{x)] 2 - A [ x[x + jx + V'{x)}) 

J Jr(0)=r o \ 2 4L) J Q J Q J 

(48) 

After some simple manipulations, one can recast this into 

7? ro (A) - e^ 1 -") J driP^Cn.Tlro.O) x cxp f-^L(A + l - >Jl + 2\)(E 1 - E )^j (49) 



where JSj = ^i 2 + (we recall A = 2DX/^,T] = V 1 + 2A). This formula is use- 

ful, since the propagator P V1 goes exponentially fast to the equilibrium value P eq (ji) = 
\fyr\j {2-kD)U~ 1 exp(—-yr]Ei/D) where U a = J dxo exp(—^faV(xo)/D) is the configurational 
partition function. Consequently, for each fixed value of ro, one has the true equivalence 



^ ro(A) ~ \[^D U ^ eUl ~ V) J driex V (-^ 5 [(\ + l + V )E 1 -(\ + 1 - V )E ]) (50) 

and one sees that there is no problem of convergence in the integral for any real value of A 
such that A > —1/2 (for A < —1/2 there is always a cut due to the presence of a square root 
in rj). Thus, one extracts exactly the leading exponential term from the preceding formula 

as 

log7f ro (A)~|(l-»7) (51) 

and shows in the same time that the large deviation function of 7r ro (e) is 

f ro (e)=-^(e-lfxe(e) (52) 

irrespective of the precise form of the potential. 

Let us now look at 7r(e): its characteristic function can be written 

t?(A) = \\ 2d\ x J dr o dv 1 P 7ri (r 1 ,r\r o ,0) x exp (-X[( A + 1 - r] )E 1 - (A - 1 - V )E ] 

(53) 

A priori, one must take care of the fact that the equivalence P(ri, r|ro, 0) ~ P e q(^i) is 
not reached uniformly with respect to r , as already noticed. But, if one inspects the free 
Brownian case, for which the exact result is computed, this replacement is finally equivalent 
to neglect exponentially small corrections (terms like tanh?77 replaced by 1 for instance); 
this is precisely this approximation which leads to formulae ( ^7|j38|) from (|35| , |36| ): the only 
limitation is that the resulting formulas are not uniformly valid in the e space. But the 
associated large deviation function is unaffected by these corrections. 
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We can assume that this scenario is still correct in the general case, an assumption which 
is very reasonable indeed. Thus, up to exponentially vanishing factors, one has 



^ (A) ~ Si X 7^77" / dr ° dri 6XP (r2D [{X + 1 + v)El ~ (A ~ 1 ~ V)E ° ] ) (M) 

= / 47? " x el ^ x U C^ + rm U (n + i-D/2 (55) 
V + - A 2 U!U V { ' 

and it is easily seen that, again, the ro integral diverges when A approaches = 4, pointing 
out the probable beginning of a real cut, already encountered when V = 0. Noticing that 
the leading exponential term is exp(^-(l — 77)), also insensitive to the presence of a pinning 
potential, one deduces that also in this case the large deviation function is given by (f^). 
Of course, fully mathematical precision is not given here, but we think that convincing 
arguments are nevertheless given in favour of our result. 

Concerning the prefactor of the large deviation function, it is interesting to mention that 
it keeps in general a dependence on the form of the pinning potential, through the functions 
U. In fact, formulas like ( p7| , p8| ) can themselves be simplified; for instance 

^»~ U--i)(.- + D» )' (56) 

(we did not propose this equivalence previously, for it is not correct in the vicinity of £ = 1/3, 
unlike formula ( |3^ ) . . . ). For the general case with a potential V, this "supersimplification" 
gives 

T(?\ „ ( ^ \ 5 v C fe- 1 + l) a /4^(3-g-i)(l+g-i)/4 ft5r > 

J{£) \n(3e — l)(e + 1) 3 / X U^Z (5?j 

Thus, the pdf associated with different potentials do not exactly coincide at large r, except 
at the value £ = 1. But, at the level of the large deviation function, the universality of /(e) is 
reached. The combination of these two points explains probably the remarkable coincidence 
of the different curves on the figure |?|. 



4 Conclusion 

In this paper, we considered the Langevin equation as a dynamical evolution of a simple 
dissipativc system driven by an external forcing, and computed the probability density func- 
tion of the time-averaged injected power in the permanent and non permanent regimes. We 
showed that the associated large deviation functions are different, in particular a negative 
tail exists only in the permanent regime. We explained the origins of this discrepancy and 
highlighted the role of the rare but very energetic initial conditions; we showed also that the 
system considered here does not verify the so-called Fluctuation Relation /(e) — /(— e) = je, 
even in the vicinity of e = 0, indicating thence that any attempt to enlarge careless the 
applicability of the conclusions of the Fluctuation Theorem to dissipative systems is vain. 

We considered thereafter Langevin equations with pinning potential, and showed that 
the associated large deviation functions are completely insensitive to the potential (but not 
the pdf itself): this result appears to be a good indication that large deviation functions 
could be an appropriate tool to characterize well general properties of systems beyond some 
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peculiar irrelevant details which faded away through the process of averaging. In our case, 
the function / tells us something global associated with the energy transfer throughout the 
Brownian particle, irrespective to the precise dynamics of each geometrical configuration. 

A natural extension will be to consider stochastic coloured dynamics, i.e. noises with 
time correlation. We want to test the "solidity" of / with respect to correlations arising 
in the external forcing. Moreover, we can address the question of the generality of the 
observed non-analyticity of the large deviation function in systems in permanent regime 
where the velocity field is unbounded: from our analysis of the Langevin system, we think 
that such anomalies could be widely encountered in dissipative far-from-equilibrium systems. 
An interesting perspective is besides to consider extremely correlated noises, in order to check 
in this context the ideas exposed in |l5[ |l6| , where some conjectures are made on the limiting 
behaviours of the pdf of global variables in highly correlated systems. 
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